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I.  Introduction 

)d 

Detectors  and  estimators  which  utilize  signals  in  the 

high  frequency  band  (HF)  as  input  data  usually  suffer  from 

performance  deterioration  as  a  result  of  the  presence  of 

high  frequency  interference.  This  interference  can  be 

* 

categorized  as  either  intelligent  or  non-intelligent. 

The  non-intelligent  interference  occurs  mainly  as  the 
result  of  channel  noise,  which  is  primarily  inpulsive  in 
nature  in  the  HF  band  (atmospheric  and  man-made  noise) 
and  should  be  accounted  for  during  the  receiver  design 
stage.  The  intelligent  interference  is  developed  by 
another  user  either  involuntarily  when  he  utilizes  the  same 
channel,  or  deliberately  to  degrade  the  receiver  performance 
The  second  type  of  interference  in  general  has  a  more 
deleterious  affect  on  the  receiver  performance  than  the 
first  and  should  be  taken  into  account  in  the  process  of 
signal  design. 

The  classical  approach  for  solving  this  problem  is 
either  to  design  a  receiver  and  to  choose  the  signal  which 
is  optimal  for  a  channel  model  developed  apriori  from 
channel  measurements,  or  to  utilize  an  adaptive  technique 
which  works  for  a  finite-dimensional  class  of  channel  models 
A  dangerous  property  of  the  first  approach  is  that  the 
performance  of  these  optimal  receivers  usually  deteriorates 
very  badly  if  the  actual  data  did  not  follow  the  assumed 


channel  model  exactly  [1]  and  [2] .  On  the  other  hand, 
continuous  changes  in  channel  models  in  HF  impulsive  noise 
channels  is  a  known  fact  for  communication  engineers,  at 
least  for  the  non-intelligent  part  of  the  noise  [3] ,  [4]  , 
and  [5].  Moreover,  there  is  no  presently  known  procedure 
to  protect  the  receiver*  completely  from  the  presence  of 
bad  data,  especially  those  due  to  measurement  errors 
either  in  the  process  of  channel  modeling  or  while  taking 
the  observations. 

The  main  difficulty  with  the  adaptive  approaches  in 
general  is  that  the  communication  system  designer  never 
knows  a  priori  which  parameters  in  his  assumed  channel 
model  should  be  allowed  to  change  with  the  measurement  of 
the  current  data. 

As  a  result  of  the  above  observations,  the  development 
and  utilization  of  new  techniques  for  the  design  of  robust 
receivers ,  which  depend  only  on  partial  information  about 
the  channel  and  provide  good  performance  over  all  the 
infinite  dimensional  class  of  channels  which  possess  the 
properties  described  by  this  partial  information,  has 
received  considerable  attention  during  the  last  few  years 
[61,  [7],  [8],  and  [91.  This  work  deals  with  the  problem 
or  receiver  design  in  the  presence  of  non-intelligent 
interference.  In  [6]  and  [7]  the  problem  of  receiver  design 
in  the  presence  of  measurement  errors  has  been  treated. 

In  [9]  a  general  approach  for  receiver  design  when  only 


partial  information  about  the  channel  is  available  has  been 
introduced.  The  results  achieved  in  this  direction  are 
highly  encouraging,  and  invite  more  effort  and  research 
toward  the  development  of  similar  techniques  to  overcome 
the  effects  of  intelligent  interference  on  the  system 
performance.  Specifically,  it  would  be  of  interest  to  both 
the  communicator  and  the  interferor  to  know  about  the  worst 
case  interference  and  the  worst  case  performance  of  certain 
receivers  and  certain  signal  forms. 

This  information  could  be  used  by  the  communicator  to 
choose  the  best  signal  waveform  and  detection  strategy 
in  the  presence  of  the  worst  case  interference.  In  other 
words,  to  design  a  system  which  maximizes,  over  all 
combinations  of  signal  waveforms  and  detection  strategies, 
the  worst  case  performance.  The  interferor  can  use  this 
information  to  prepare  an  interference  strategy  which  gives 
maximum  degradation  of  the  system  performance  without 
apriori  knowledge  of  the  signal  waveform  at  each  point  of 
time  and  within  his  constraints,  such  as  transmitted  power 
limits. 

As  a  compatible  effort  in  this  line  of  research  the 
test-bed  package  "Detector  Performance  Evaluation"  (DPE)  was 
designed  to  evaluate  the  performance  of  a  wide  class  of 
detectors  for  different  signal  waveforms  against  arbitrary 
noise  distributions.  A  detailed  description  of  the  contents 
of  this  package  is  given  below. 


-4- 


II .  Package  Description 

The  test-bed  package  DPE  in  a  fortran  oriented  simulation 
package  which  consists  of  three  programs  and  two  supporting 
libraries.  These  three  programs  are  DPEP.FOR,  DPPRE . FOR, 
and  DPPOST . FOR.  The  two  libraries  are  called  DPED  and  DPESN. 
The  main  program  in  the  .package  is  DPEP.FOR  which  runs  with 
the  support  of  the  above  two  libraries.  The  other  two 
programs  run  independently.  No  other  libraries  or  subroutines 
are  required  other  than  the  standard  fortran  library.  The 
first  program  to  be  executed  should  be  DPPRE. FOR.  The 
main  function  of  this  program  is  to  prepare  an  input  file 
for  the  main  program  DPEP.FOR.  The  name  of  this  file  is 
DP IN . DATA  and  it  contains  all  required  information  about 
the  detector  to  be  tested,  the  signal  waveform  and  the  noise 
distribution. Utilization  of  the  program  DPPRE  facilitates 
running  the  program  DPEP  off-line.  To  limit  the  storage 
requirements,  the  output  of  DPEP  will  be  a  binary  file 
which  can  be  translated  into  a  formated  file  using  the 
program  DPPOST.  The  library  DPED  contains  10  detector 
functions,  five  of  which  are  already  specified  and  five  of 
which  are  left  to  the  user  to  specify.  These  detector 
functions  are  called  dl,d2,...,  dlO.  The  library  DPESN 
contains  both  signal  waveforms  and  noise  distribution 
functions.  It  can  accept  up  to  5  different  waveforms  called 
sl,s2,...,  s5  and  10  different  noise  distributions  called 
znl,zn2,...,  znlO.  Both  these  libraries  are  extendable. 

A  detailed  description  of  these  programs  and  libraries  is  to 
be  given  below. 


II. 1  DPPRE. FOR 


The  main  function  of  this  program  is  to  prepare  an  input 
data  file  for  the  main  program  DPEP.FOR.  The  inputs  to  the 
program  are  the  answers  to  a  set  of  questions  which  appear 
on  the  terminal.  The  output  of  the  program  is  a  file  called 
DP  IN .  DATA  of  format  (5i^>,  5fl0.6).  This  file  contains  the 
previous  answers  in  the  sequence  required  by  the  program 
DPEP.FOR  which  utilizes  this  information  to  define  the 
detector  to  be  tested,  the  signal  waveform  and  the  noise 
distribution  function  and  parameters.  Beside  each  question 
the  format  of  the  answer  will  appear  to  the  user. 

The  first  question  to  appear  is  of  the  form  "detector 
no.  =  ?".  The  answer  to  this  question  is  an  integer 
between  1  and  10  and  defines  which  detector  function  in 
the  library  DPED  is  to  be  tested.  For  example  if  the  answer 
is  3,  then  the  detector  function  to  be  tested  is  d3  which 
is  the  "Limiter-Correlator  Detector". 

The  next  two  questions  define  the  noise  distribution 
function.  The  first  one  will  be  of  the  form  "noise  =  ?". 

The  answer  to  this  question  defines  the  noise  distribution 
l  nction.  For  example,  if  it  is  1  then  the  noise  distribution 
function  will  be  znl  which  is  the  normal  distribtion  in 
the  library  DPESN.  The  answer  to  this  question  should  be  an 
integer  between  1  and  10.  The  answer  to  the  second  question 
which  appears  in  the  form  "scale  =  ?"  affects  the  program 
DPEP  in  different  ways  depending  on  the  answer  to  the  previous 
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question.  If  the  answer  to  the  previous  question  was  1  or  2 
then  the  answer  to  this  question  determines  the  noise 
standard  deviation.  If  the  previous  answer  was  3  then  it 
determines  the  standard  deviation  of  log(x).  If  the  previous 
answer  was  4  or  5  then  this  number  defines  the  standard 
deviation  of  the  contaminating  distribution;  the  standard 
deviation  of  the  contaminated  distribution  is  always  one. 

The  next  two  questions  to  appear  are  signal  related 
questions.  The  first  is  about  the  signal  amplitude  and 
has  the  form  "sig.  amp.  =  ?" .  The  second  is  about  the 
signal  waveform  and  has  the  form  "sig  =  ?".  The  answer 
to  this  question  should  be  an  integer  between  1  and  5,  and 
defines  which  signal  waveform  from  those  in  the  library 
DPESN  should  be  used.  This  waveform  is  also  utilized  as 
the  reference  signal  for  correlation  type  detectors. 

The  next  two  questions  define  both  the  sample  size  and 
the  number  of  runs  to  be  averaged.  They  are  of  the  form 
"sample  size  =?  "  and  "no.  of  runs=1000x?"  respectively. 

The  answer  to  the  first  one  should  be  a  positive  integer 
less  than  or  equal  to  30.  The  second  can  take  any  integer 
value  acceptable  by  the  machine. 

The  next  set  of  questions  has  the  form  "Upper  Threshold=?" 
and  "Lower  Threshold=?"  and  determines  the  threshold (s)  of 
the  test  and  the  type  of  test  to  be  performed.  If  the  upper 
threshold  (tu)  is  equal  to  the  lower  threshold  (tl) ,  then 
it  is  a  one  threshold  test,  and  the  program  DPEP  determines 
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the  probability  that  the  test  statistic  is  on  either  side  of 
the  threshold  as  is  the  case  in  the  minimum  probability  of 
error  criterion.  If  tu  is  not  equal  to  tl ,  the  distance 
between  tl  and  tu  will  be  divided  into  500  equal  intervals, 
and  the  results  of  DPEP  will  be  in  the  form  of  a  histogram. 

If  the  answer  to  question  1  in  this  program  was  4  or 
5  then  tl  and  tu  should  be  between  +8.0. 

The  last  question  here  is  the  form  "a=?"  and  it 
appears  only  if  the  answer  to  question  1  was  5.  The  answer 
here  defines  the  break  point  for  the  nonlinearity  utilized 
by  detector  d5.  More  details  for  this  number  will  be  given 
later. 

II. 2  DPEP . FOR 

This  is  the  main  program  in  the  test-bed  package  DPE. 

It  consists  of  four  routines,  the  main  routine  and  three 
subroutines  d,  c  and  e.  It  runs  with  the  support  of  the 
two  libraries  DPED  and  DPESN.  The  input  to  this  program 
is  the  file  DPIN.DATA  which  is  the  output  of  the  program 
DPPRE .  The  output  file  is  an  on-record  binary  file  of 
length  (2012)  which  contains  the  number  of  observation 
groups,  the  upper  and  lower  thresholds,  and  the  number  of 
runs  in  which  the  test  statistic  falls  in  some  interval. 

The  main  routine  and  the  two  subroutines  d  and  c 
define,  respectively,  for  the  subroutine  e  which  detector, 
noise  distribution,  and  waveform  combination  is  to  be 
tested.  The  actual  performance  eavaluation  test  occurs 
in  subroutine  e. 


This  subroutine  utilizes  the  function  call  name  delivered 
to  it  by  the  three  other  routines  and  the  rest  of  the 
information  in  the  input  file  to  calculate  the  value  of  the 
test  statistic  under  the  assumption  of  independent  data. 

If  the  test  isaqne  threshold  test  it  compares  the  test 
statistic  with  the  threshold.  If  it  is  a  two  threshold 
test,  it  divides  the  distance  between  the  two  thresholds 
into  500  equal  intervals  and  decides  in  which  interval  the 
test  statistic  falls,  and  adds  1  to  the  counter  for  this 
interval.  Every  one  thousand  runs,  this  subroutine 
updates  the  output  file. 

II. 3  DPPOST. FOR 

This  program  translates  the  binary  file  DPEP . DATA , 
which  is  the  output  of  the  program  DPEP,  into  a  formatted 
file  called  DPOUT.DATA.  The  first  record  in  this  file 
contains  the  number  of  runs.  The  second  contains  two 
columns ;  the  first  contains  the  threshold  and  the  second 
contains  the  probability  that  the  test  statistic  exceeds 
this  threshold  (the  probability  of  false  alarm  or  detection  depending 
on  whether  the  signal  is  present  or  not) .  A  sample  of  this 
file  is  shown  on  the  next  page. 
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nunber  of  runs  ■ 
Threshold  <  T  > 
1  . 7  50  0  C 
1 .76000 
•  1 .77000 
1 .78000 
1.79000 
1 .80000 
1 .81000 
1 .82000 
1 .83000 
1 .84000 
1  .85000 
1 .86000 
1 .87000 
1 .88000 
1 .89000 
1 .90000 
1 .91000 
1 .92000 
1 .93000 
1 .94000 
1 .93000 
1 .96000 
1 .97000 
1 .98000 
1 .99000 
2.00000 
2.01000 
2.02000 
2.03000 
2.04000 
2.05000 
2.06000 
2.07000 
2 .08000 
2.09000 
2.10000 
2.11000 


2000 .0 

p<  x>T > 
0.78799999 
0.77550000 
0.76400000 
0.75349998 
0.74349999 
0.73600000 
0.72600001 
0.71799999 
0.70749998 
0.69900000 
0.69250000 
.  0.68049997 

*  0.67150003 

0.66149998 
0.65050000 
0 .63700002 
0 .62449998 
0.61350000 
0.59700000 
0 .58550000 
0.37349998 
0.35849999 
0 .34449999 
0.53330002 
0.51349998 
0.50000000 
0 .48600000 
0.46950001 
0.45899999 
0.44450000 
0 .43000001 
0.41400000 
0.40200001 
0.39050001 
0.38000000 
0 .37349999 
0.36300001 
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ll .  4  The  Library  DPED 

This  libary  contains  up  to  10  detector  functions. 
Five  are  specified  and  are  briefly  described  below. 
Throughout  this  section  we  shall  consider  as  the  ith 
observation,  n  as  the  total  number  of  observations, 
as  the  ith  signal  samp*le,  and  Tn  as  the  test  statistic. 


II. 4.1  function  dl  (Linear  detector) 

n 


n 


1  11 
=  -  l 

n  i=l 


X. 

l 


II. 4. 2  function  d2  (correlator  detector) 

n 


n 


=  i  l 

n  .L. 


i=l 


Si  Xi 


II. 4. 3  function  d3  (Limiter  -  correlator  detector)  [10] 


l  n 

p  =  -  y  s.  Kx.) 

n  n  1  1 


where 


1(X.)  = 


-1.14 

X. 

l 

1.14 


X^-1.14 

|X.|<1.14 

X . >1 . 14 
1— 


II.  4. 4  function  d4  (M-detector  for  conteiminated  normal  class)  [9  ] 


T„  is  such  that 

N 

N 


Si  1(VTNSi>  -  ° 


1 ( . )  is  the  same  as  for  the  previous  detector 
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II. 4. 5  function  d5  (M-detector  for  p-point  class)  [9  ] 


T„  is  such  that 
N 


J,  Si  1(X.-V.)  =  0 


If  the  class  is  defined  by 


F  =  {f  :  /  f (x) dx  <  1/2} 


then 


-tan  (ca) 
tan (ct) 
tan  (cot) 


t<-a 


1 1  <a 


1.718a 


Notice  that  the  parameter  a  is  the  same  as  a  in  the 
program  DPPRE. 


II. 5  The  Library  DPESN 


This  library  consists  of  5  signal  waveforms  and  10 
noise  distribution  functions.  Of  the  five  signals  only  one 
is  specified  as  a  constant  signal.  Of  the  ten  noise 
distribution  functions  the  following  five  are  specified. 

The  parameter  an  in  all  these  functions  is  the  scale  factor 
in  the  program  DPPRE. 


.  ,  illfc  (  _ 
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II. 5.1 

function 

znl 

(normal  distribution) 

f  (x) 

1 

(x2)  1 

a  /Tn  exp 

n 

L  2  a  2  J 

n 

II. 5. 2 

function 

zn2 

(double  exponential) 

f  (x) 

= 

A 

2a- 

n  L 

M  I 

an  J 

II. 5. 3 

function 

zn3 

(lognormal  component) 

x  =  y 

cos 

0 

f  (y) 

1  I 

Qvr>  i 

log2 (x) 

]  yi° 

_ _  ***■?  i 

xa  ZTn 

n  1 

2a  2 

- 

f  (0) 

=  . 

1 

2tt 

|e|<n 

II. 5. 4 

function 

zn4 

(normal  contaminated  by 

normal) 

f(x) 

-x2 

9  _1  e 

+  1  1 

-  X2 

2an^ 

/Tn 

a  /Tn 
n 

G 

II. 5. 5 

function 

zn5 

(normal  contaminated  by 

double  exponential) 

f  (x) 

= 

2 

-x 

.9  —  e  2 

/Tn 

+  -1  IS"  6 

n 

-M 

a 

n 
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III. Updating  and  Extending  the  Package 

1.  The  maximum  sample  size  allowed  in  this  package  is  30. 

To  increase  this  limit,  the  user  should  change  the 
dimension  of  the  arrays  xd  and  ref  in  instruction 
number  2  in  the  subroutine  of  the  program  DPEP. 

2.  To  add  a  new  detector  function  —  within  the  limit  of 

a 

ten  detectors  —  the  user  has  only  to  replace  one  of 
the  detector  functions  in  library  DPED  by  his  new 
function.  If  we  want  to  replace  the  detector  function 
d^  (where  i  is  an  integer  between  1  and  10) ,  the  call 
name  for  the  new  detector  function  should  be 
d^ (xd,n,y , z) ,  where  xd  is  an  array  for  the  samples,  n 
is  the  sample  size,  y  is  an  array  for  the  reference 
signal  samples,  and  z  is  an  optional  variable. 

3.  Similar  methods  would  be  used  to  update  the  signal 
waveforms  and  the  noise  functions  in  the  library  DPESN . 

4.  To  extend  the  number  of  detector  functions  beyond  the 
limit  of  10,  the  user  must  change  both  library  DPED  and 
the  program  DPEP.  For  example  if  the  user  wants  to  add 
detector  number  11,  then  a  detector  function  dll(x,n,y,z) 
has  to  be  added  to  the  library  DPED.  Then  a  statement 
of  the  form,  if(i  .  eq.  11)  call  dll (idata,dat,dll) 

has  to  be  added  just  before  the  "END"  statement  in 
the  main  routine  of  the  program  DPEP.  The  external 
statement  in  this  routine  has  to  be  changed  to  contain 


dll. 
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5.  To  extend  the  number  of  signal  waveforms  beyond  the  limit 
of  5,  the  user  must  add  his  new  waveform  to  the  library 
DESN.  Assuming  it  is  waveform  number  6  then  the  call 
name  for  the  new  function  would  be  S6(j),  where  j  is 
the  order  of  the  signal  sample.  The  next  step  is  to 
add  a  statement  of  the  form 

if  (is  .eg.  6)  call  e (idata,data,x ,y ,S6) 
just  before  "RETURN"  in  the  subroutine  "c"  of  the 
program  DPDP.  The  second  statement  in  this  subroutine 
must  be  changed  also  by  adding  S6  to  the  set  of  external 
functions;  i.e.,  this  statement  should  be  changed  from 
external  SI ,S2 ,S3 ,S4 ,S5 

to 

external  SI ,S2 ,S3 ,S4 ,S5 ,S6  . 
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dppre.for  -  A  progrtn  to  prepare  input  data  file  for 
the  na i n  program  (dpep.for)  {Detector 
Performance  Evaluation  Progran). 


Inputs  -  A  set  of  nunbers  to  be  delivered  through  (tty) 
as  answers  to  questions.  The  fornat  of  each 
number  appears  next  to  the  question.  This  set 
of  nunbers  determines  the  detector  to  be  tested* 
the  noise  distribution  function*  the  signal* 
and  the  thresholds. 


Outputs  -  A  file  which  contains  all  the  above  infornation 
and  serves  as  an  input  for  dpep.for.  File 
fornat  is  < 51 5 » 5f 1 0 . 6 > .  The  name  of  this  file 
is  <  dp i n  .  da  ta  ) . 


0001 

0002 


dinension  i  data<  5  >*  data<  5  > 
call  setfil  ( 1 *' dpi n . data' ) 


This  section  deternines  the  detector  to  be  tested.  The 
answer  to  this  question  is  an  integer  between  I  and  10 
detector  no.  detector's  nano 


0003 

0004 

0005 

0004 


5  . 

£  -  10 


Linear  detector. 

Correlator  detector. 

L i n l ter-Corr el ator  detector. 
M-detector  for  contaminated 
normal  class. 

H-detector  for  P-point  class. 
Opt i ona l . 


write  (4*1) 
fornat  ('detector  no. 
read  (5*2)  idata(  1) 
fornat  (i5> 


■?  ( i3>' > 


Tihs  section  deternines  the  noise  dltrlbution 
function.  The  answer  to  the  next  question  should  be 
an  integer  between  1  and  10  according  to  the  following 
table. 


Distribution  no 
1 
2 
3 


0 i str i but  ion 
Nornal . 

Double-Exponential. 

A  component  of  a  Log¬ 
normal  distribution. 
Contan i nated-Nornal * 
the  contaminating  dist 
is  also  nornal  with 
optional  variance.  The 
contamination  ratio  is 
lOt. 

Contaminated-Normal  as 
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c  in  4  exept  that  the 

c  contan I  not i ng  distribution 

c  is  a  Doub le-Exponenti a l . 

c 

0007  writ*  (8.8) 

0008  8  .  fornat  ('non**?  (15)'  > 

0009  read  (3,2)  idata<2> 

c 
c 

c  The  answer  to  the  next  question  should  be  a  reat 
c  positive  nunber.  If  the  noise  is  of  type  I  or  2.  then 
c  this  nunber  deternines  the  standard  deviation.  If  the 
c  noise  is  of  type  3.  then  it  deternines  the  standard 
c  deviation  of  log(x).  If  the  noise  is  either  of  type 
c  4  or  5.  then  this  nunber  deternines  the  standard 
c  deviation  of  the  contan inati ng  di str ibut i on . 
c 

0010  write  (8.93) 

0011  93  fornat  ('scale*?  <fl0.8)') 

0012  read  (5,94)  data( 1 > 

c 
c 

c  The  answer  here  should  be  a  real  nunber  which  gives 
c  the  signal  anplitude. 
c 

0013  write  (8, 10) 

0014  10  fornat  ('sig.  anp.*?  (fl0.8>'> 

0015  read  (3/94)  data<2> 

c 
c 

c  This  is  to  define  the  signal  wave  forn.  If  the 
c  answer  is  <  2.  the  signal  will  be  considerd  as  a 
c  constant.  The  sane  signal  waveforn  will  serve  as 
c  a  reference  signal  for  correlator  and  H  detectors, 
c 

0018  write  (8.3) 

0017  3  fornat  C'sig-7  (15)') 

0018  read  (3.2)  idata(3) 

0019  if  (idata(3)  .le.  1>  idata(3)*l 

c 
c 

c  This  is  to  deternine  the  sanple  size.  It  should 

c  be  an  integer  between  1  and  30. 
c 

0021  write  (8/90) 

0022  90  fornat  ('sanple  size"?  (15)') 

0023  read  (3*2)  idata<4> 

c 

c  This  is  to  deternine  the  nunber  of  thousands 
c  of  runs, 
c 

0024  write  (8*92) 

0025  92  fornat  ('no.  of  runs  *1000e?  (IS)') 

0028  read  (3.2)  idata(S) 

c 

c 

c  This  section  deternines  the  threshold  (s)  of  the 
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c  test  and  the  type  of  test  to  be  perforned.  If  the 
c  upper  threshold  <tu>  is  equal  to  the  lower  threshold  <tl>/ 
c  then  >t  is  a  one  threshold  test/  and  the  progran  (dpep> 
c  determines  the  probability  that  the  test  statistic  is  in 
c  either  side  of  the  threshold  as  is  the  case  in  the 
c  n ini  nun  probability  of  error  criterion.  If  tu  is  not 
c  equal  to  tl<  the  distance  between  tl  and  tu  will  be 
c  divided  into  300  equal  intervals/  and  the  results  of 
c  (dpep)  will  be  in  the  forn  of  a  histogram, 
c 

0027  write  <6/93) 

0028  93  format  ('Upper  Threshold  *?  <fl0.6>'> 

0029  read  (5/94)  data<3> 

0030  94  format  (fl0.6) 

0031  write  (6/97) 

0032  9 7  format  ('Lower  Threshold  *?  <fl0.6>') 

0033  read  (3/94)  data(4> 

0034  tref*0.Q 

0033  if  (idata(l)  .  ne  .  5)  go  to  99 

c 
c 

c  This  number  determines  the  break  point  of  the 
c  nonlinearity  in  detector  no.  3.  This  is  the  "a*  parameter 
c  in  the  class  of  8-point  distributions. 


c 


0037 

write  (6/98) 

0038 

98 

format  <'a«?  (flO.C)'  ) 

0039 

read  (5/94)  data<5> 

0040 

99 

write  (1/100)  (Idata(i)  /i»l/ 

0041 

write  (  1/200)  <data( 1  )  /  i »1 / 5  ) 

0042 

100 

format  (SiS) 

0043 

200 

format  <5fl 1 .3  > 

0044 

end  file  1 

0045 

end 
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DPEP.F0R-*  0  progran  to  evaluate  the  perfornance 
of  detectors  in  the  presence  of  different 
signals  end  against  optional  noise  densities. 

This  progran  includes  in  addition  to  the  nain 
routine  the  three  subroutines  d(.)«c<.)  and  e<  .  )* 
and  utilizes  different  functions  fron  the 
tuo  libraries  <dped)  and  (dpesn).  The  first 
tibrory  contains  all  detectors  to  be  tested.  The 
second  library  contains  different  signals 
and  noise  generators. 


Input>  File  <dpin.data>  which  is  the  output  of  the 
data  preparation  progran  (dppre.for).  File 
f ornct  is  <  5 iS,3f 10 .6  > . 


Output:  o  binary  file  called  <dpep.data>.  This  file  consists 
of  one  record  of  length  2012  and  contains  the  nunber 
of  observation  groups  utilized*  the  upper  and  lower 
thresholds  and  the  probability  of  false  alarn  or 
detection  depending  on  whether  the  signal  is  present 
or  not.  This  file  can  be  transforned  to  a  fornated 
file  using  the  progran  < dppost . f or >. 


Wain  routine: 

dl #d2»  . . . #dl0  is  a  set  of  external  functions 

fron  the  library  <dped)  and  defines  the  detector 

to  be  tested  according  to  the  following  list: 

dl  linear  detector. 

d2  correlator  detector. 

d3  l i n  iter-corre lator  detector 

d4  H-detector  for  contaninated  nornal  class 

d5  n-detector  for  p-polnt  class 

d6~dlQ  optional 

Depending  on  the  value  of  (■idataf!)#  the  nain 
routine  calls  the  subroutine  d< i data < data • di  ) . 


0001 

t  x  t  #  r  f)  o  I  dl  /  d2/d3/d4*  dS  /  d 6  *  d?  •  d 8  /  d 9 »  dlO 

0002 

d i nens i on 

i  data(  5  ># 

data( 5) 

0003 

call 

setf  i 

l  <2  » *  dpi 

n. data' > 

0004 

read 

(2,99)  <  i  da ta< 

i  )i  i » 1 •  5  > 

0005 

read 

(2.97)  <  data ( i 

>«  i*l ,5  > 

0000 

97 

f ornat  < Sf 1 1  .5 > 

0007 

91? 

fornct  <5» 

5  ) 

0008 

i - i da ta( 1 > 

0009 

if  <  i 

♦9 

1  >  call 

d(  i  data  /  data  *  dl  > 

0011 

iff. 

.eg 

2  >  call 

d<  i data#  data •  d2  > 

0013 

if  <  i 

.eg 

3)  call 

d< i data  #  data  <  d3  > 

0015 

if  <  i 

.  eg. 

4  >  call 

d< i data  # data • d4  > 

0017 

If  C  i 

•  *9  • 

5  >  call 

d< i data ,  data •  d5  > 

0019 

1  f  <  l 

eg. 

6  >  call 

d< i data  # data •  dt > 

j 
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0021 

if  (  »  . 

«q 

?  > 

ca  l  1 

d(  i  da  t-a  /  data  /  d?  > 

0023 

if  (  i  . 

*q 

8  > 

call 

d(  i data  >  data  >  d8  ) 

0025 

If  <  i  . 

*q  • 

9  ) 

ca  1 1 

d< i data • data / d9  ) 

0027 

If  <  1  . 

•q 

10> 

call 

dC  i data / data  <  d  10  ) 

0029 

end 
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c 

0001  subroutine  d(  i  do  t  a ,  da  t  a ,  x  > 


Subroutine  d( i da ta / da ta • x > > 

idato  and  dote  ore  too  arrays  fron  the  no i n 
routine. 

x  is  a  dunny  argunent  to  be  replaced  by  the 

coll  none  of  the  detector  under  test. 

xn 1 /  zn2 i . . . •  zn 10  is  o  set  of  external  functions 

fron  the  library  (dpesn)  and  defines  the  noise 

distribution  to*be  utilized  according  to  the 

following  list. 

znl  nornol  distribution 

zn2  douple  exponential 

zn3  a  conponent  of  a  log-nornal  distribution 
zn4  centaninated  nornol/  the  contan i not  i  ng 
distribution  is  also  nornal  with 
optional  variance.  The  contani not  ion 
ratio  is  102. 

zn9  contani nated-norna  l  as  zn4  except  the 
contani nat < ng  distribution  is  a 
doubt  e-exponent  i  a  l  . 
zn6-znl0  optional 

Depending  on  the  value  of  j«idata(2>/  this 
subroutine  calls  the  subroutine  c< i da  to • da ta • x / zn< j  >  ) 


external  znl<zn2/zn3/zn4/zn9/zn6.zn7/zn8/zn9/znl0 
dinension  I data< 1  )/ da ta< 1  ) 
i n«i data(  2) 

if  < in  .eg.  1>  call  c(  i  data  /  data / x / zn  1  > 
if  (in  .eg.  2>  call  cC i data # data / x / zn2 ) 
if  < >n  .eg.  3)  call  c<  i data / dot  a . x • zn3 ) 

if  (in  .eg.  4)  call  c(  idata/data#x/zn4) 

if  (in  .eg.  9)  call  c(  idata / data / x / zn9  ) 

if  (in  .eg.  6)  call  c( i data • data / x # znt > 
if  (in  .eg.  7>  calt  c(idata/dota.x.zn7) 
if  (in  .eg.  8)  call  c(  i  data / data / x / zn8 ) 

if  (in  .eg.  9>  call  c<  i  data / data / x / zn9 > 

if  (in  .eg.  10>  call  c< i data/ data/ x/  znlO) 
return 


0002 
0003 
0004 
0009 
0007 
0009 
001  1 
0013 
0019 
0017 
0019 
0021 
0023 
0029 
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0001 


subroutine  c<  i  data*  data*  x*  y  ) 

Subroutine  c  < I data* data* x*y > 

idata*  end  data  are  two  arrays  fron  routine 
nain  through  routine  d(. *.#.). 
x  is  a  during  argunent  to  be  replaced  by  the 
call  nane  for  the  detector  under  test, 
y  is  a  dunny  argunent  to  replaced  by  the 
call  nane  for  the  noise  distribution  function 
fron  subroutine  <} . 

sl«s2*...*s5  is  a  set  of  external  functions 
fron  library  (dpesn)  and  defines  the  signal  vave 
forn.  It  also  defines  the  reference  signal 
for  sone  detectors.  si  is  a  constant  signal  and 
others  are  optional. 


Depending  on  the  value  of  k*idata<3>/  this 
subroutine  calls  the  subroutine  e< I  da ta * da ta • x * y • sk > . 


0002 

externa  l 

si  * 

s2  # 

s3«s4#  s5 

0003 

d  i  net's  i  on  idata(l)/ 

data(  1  > 

0004 

i s» »  da  t  a<  3  > 

0005 

i f  ( is  . 

•  q- 

1  > 

call 

e(  idata«data/X*y/Sl > 

0007 

i f  < is  . 

•q- 

2) 

call 

e< idata*data*x*y<s2> 

0009 

i  f  <  i  s  . 

•  q- 

5  > 

call 

e< idata*data*x/y»s5) 

001  1 

i  f  C  i  S 

•q- 

4  ) 

call 

e< idata«data/X>y/s4) 

0013 

if  (is  . 

•  q- 

3) 

call 

e( idata*data>x/y*s3) 

0015 

return 

0016 

end 
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0001 


subroutine  e< l date, data* x, y«z > 

Subroutine  e  < I  date • dot® •  x # y , z  ) » 

This  is  the  n«in  routine  in  the  test 
package  where  actual  detector  test  is  perforned. 

idata*data»x  and  y  are  as  described  in  the  subroutine 
c  . 

z  is  a  dunny  argunent  to  be  replaced  by  the  call 
none  for  the  sighal  wave  forn  fron  subroutine  c. 

redata(l)/  defines  the  standard  deviation  of  the 
noise  distribution  if  it  is  znl  or  zn2  and  defines 
the  standard  deviation  of  log(x>  if  the  distribution 
is  zn3.  It  defines  the  standard  deviation  of 
the  contani nat  ing  distribution  in  case  of  zn4 
and  znS. 

t*data(2)  is  the  signal  anplitude. 

tu«data(3>  and  tl«data<4>  defines  the 

upper  and  lower  thresholds  respectively.  If  tu»tl 

the  test  to  be  perforned  is  a  one  threshold  test. 

If  tu  is  not  equal  to  tl  the  output  of 
the  test  is  a  histogran  with  the  distance 
between  tl  and  tu  divided  into  500  equal 
i nterva  Is  . 

tref *data< 5 > #  is  utilized  only  with  d3  to 
define  the  break  points  of  the  non l i near i ty . 


n*idata(4)  is  the  sanple  size. 


nl*tdata<5>  Is  the  nunber  of  thousands  of  tines 
the  test  to  be  perforned. 


0002 

c 

dinension  ldata<  1  >/  datat  1  >#ox(30><tr<500>/ref<30> 

0003 

call  setfll  < 1 >  ' dpep . data' ) 

0004 

define  file  1  < 1 / 2012/ u, kx ) 

0005 

i *  ldata< 1 > 

0006 

r*da tat  1 > 

0007 

t«data< 2  > 

0000 

tu«data( 3 > 

0009 

1 1 «data< 4  ) 

0010 

tref >data(S > 

0011 

n«  idatat  4  ) 

0012 

nt«idata<5> 

0013 

99 

su*tu-t l 

0014 

do  100  n*l/n 

0015 

10C 

r  eft  n  >«zt  n  > 

001* 

do  101  n2*l*nl 

0017 

do  102  n«l»1000 

0010 

do  103  id«l«n 

0019 

103 

ox(  i  d  >*y(  r  > 
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0020 

if  <  t  .  eq .  0 .0  >  go  to  105 

0022 

do  104  id*l.n 

0023 

104 

ox<ld>"OX(id>+t*ref(id> 

0024 

105 

st»x<©x»n/ref.tref> 

0025 

if  (  sw  .eq.  0.0)  go  to  106 

002? 

if  <st  .ge.  tu >  go  to  200 

0029 

if  ( st  .It.  tl >  go  to  102 

0031 

k*int((<st-tl>/<tu-tl > )»500*1 ) 

0032 

go  to  210 

0033 

200 

k>500 

0034 

210 

tr(k  >«tr<  k)M 

0035 

go  to  102 

0034 

106 

if  (st  .gt.  tu >  go  to  107 

0030 

tr<  1  >«tr<  1  >♦  1 

0039 

go  to  102 

0040 

10? 

t  r  <  2  )■  t  r<  2  >♦  1 

0041 

102 

cont inue 

0042 

101 

write  <1'1>  n2 . t 1 » tu» < tr < )  > »  J» 

0043 

end  file  1 

0044 

return 

0045 

end 

-  . - 


1  .500  > 
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c  DPPOST .FOR 
c 

c  ft  progran  to  transforn  the  binary  file  (dpep.data) 
c  to  a  f ornated  file, 
c 
c 

c  Input: 

c  (dpep.data)  a  one  record  binary  file  of  length 
c  (2012).  This  is  the  output  of  the  progran  (dpep.for). 
c 
e 

c  Output: 

c  (< dpout . data >  a  fdrnated  file.  First  line  is  the 
c  nunber  of  runs.  The  rest  of  the  file  consists  of  two  cotunns 
c  containing  threshold  <T>  vs.  probability  of  exceeding  this 
c  threshold. 


c 

c 


0001 

dinension  t( 500 > * b< 300  ) 

0002 

call  setfil  ( 1 / * dpep. data'  > 

0003 

define  file  1  < 1  *  2012/ u< kx  ) 

0004 

call  setfil  <2 /' dpout . data '  ) 

0009 

read  (l'l)  n .  1 1 /  tu,  <  t<  i ) ,  i  « 1  *  300  > 

0004 

dn*n*1000 .0 

0007 

write  (2/1000)  dn 

0008 

1000 

fornat( 'nunber  of  runs  ■  '»fll.l> 

0009 

write  (2/1300) 

0010 

1500 

f  ornat(  '  Threshold  (T)  '  /  I3x  / '  p(  x  >T  >' ) 

0011 

C 

if  ( tu  .eq.  tl >  go  to  100 

0013 

c 

do  10  1*1/499 

0014 

J "500- i 

0013 

10 

t(  J  >*t(  J  >et(  j ♦  1  > 

0018 

dt*(  tu-tl >/300  0 

0017 

do  20  1*1/300 

0018 

t<  i  )*t(  i  >/dn 

0019 

20 

b( 1  )*t 1 ♦( i-l >»dt 

0020 

do  30  1*1/300 

0021 

30 

write  (2/2000)  b(  i  ).t(  i  > 

0022 

2000 

fornat  (f 10 . 5/ lOx #f 20 . 8  ) 

0023 

c 

go  to  200 

9 

0024 

c 

100 

t<  2)*t( 2>/dn 

0025 

wr  1  te(  2/2000  )  tl/t(2) 

0028 

200 

continue 

0027 

end  file  1 

0028 

end  file  2 

0029 

end 

UNIX  f of  tron  iv  vOl-11  source  listing 


page  001 


0001 


0002 

0003 

0004 


OPED . FOR 

A  sot  of  dotoctor  functions  to  support  the 
progran  DPEP. 


For  any  of  these  functions  xd  is  the  observation  array* 
n  is  the  sanple  size  and  ref  is  the  reference  signal 
array.  The  variable  tref  is  used  only  uith  the  detector 
function  dS  as  the  break  point  of  the  nonlinearity. 


real  function  dl<xd*n* ref * tref > 

Linear  Detector 

d i nens i on  xd< 1 > 
yd»0 .0 

do  201  il»l*n 


0005  201  yd«yd+xd<tl> 


0004 

0007 

0000 

0009 


yd«yd/n 

dl-yd 

return 

end 
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0001  real  function  d2< xd>n< ref / tref > 

c 

c  Correlator  Detector 


c 


0002 

dinar**  ion  xd(l>/ref<l) 

0003 

y  d»0 . 0 

0004 

do  21 I  ii*t«n 

0003  21  1 

yd»yd+xd< i i >*ref < i  i  > 

0004 

yd«yd/n 

0007 

d2*yd 

0000 

return 

OOOf 

end 
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0001 


0002 

0003 

0004 

0005 

0007 

0008 

0009 

0010 

0011 

0012 

0013 

0014 


real  function  d3< xd * n / re f # tre f > 
c 

c  Liniter-Corretator  Detector 
c 

dinension  xd<l>»ref(!> 
yd«0.0 

do  221  ii«l*n 

if  <«bs(xd<ii>>  .ge.  1.14)  go  to  222 
ydt*xd( i I  >eref < I »  ) 
go  to  223 

222  y d 1» 1 . 1 4*  ref < l i >•* l gn( 1 . « xd< I i  )  > 

223  ydayd+ydl 
221  continue 

d3«yd/n 

return 

end 
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0001 


0002 

0003 

0004 

0009 

0004 

0007 

0000 

0009 

0010 

0011 

0013 

0014 

0019 

0014 

0017 

0010 

0019 

0020 

0021 

0022 


real  function  d4(xd,n, ref , tref > 

c 

c  H-Detector  for  Contan mated  Nornal  Class 
c 

dmension  xd<l>,ref<l> 

.  eu*8.0 
el*-0.0 
en«0  .0 

do  231  ilsi,20 
dxt*0 . 0 
do  232  jj*l,n 
ens*en*ref< J  J  > 
dx*xd< 1 J >-ens 

if  <abs(dx)  .ge.  1.14)  dx» 1 . 14es tgn< 1 .  • dx > 
232  dxt«dxt*dx*ref < J J  > 

If  <  dxt  >  234,239,234 
234  eu«en 

go  to  237 

234  el«en 

237  en*< eu+el >/2 .0 
231  continue 

235  d4»en 
return 
end 
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0001  real  function  d3< xd*n> ref  * tref > 

c 
c 

c  N-Oetecter  for  P-Point  Class 


c 


0002 

dinension  xd<l)«rtf(l) 

0003 

*u»8 .0 

0004 

•l»*0.0 

0003 

e*»0 .0 

0004 

do  241  li«l,20 

0007 

d:<t«0 . 0 

0000 

do  242  jj*l*n 

0009 

ens«t«*ref< j  J ) 

0010 

dx"xd<  J  J  >-ens 

0011 

if  (abs<dx>  .1*.  tref>  go 

0013 

dx«0 .44 3* si gn< 1 . »  dx ) 

0014 

go  to  244 

0013 

24  3 

d  x •( dx/Cl . 7 1 8*tr e  f  ) ) 

0014 

dx«s in<  dx  >/cos(dx  > 

0017 

244 

dxt»dxt*dx*ref < j j  ) 

0010 

24  2 

continue 

0019 

if  C  dxt  >  24  3  < 244  «  247 

0020 

24  5 

eu«en 

0021 

go  to  248 

0022 

24  7 

•  l  ■M 

0023 

24  6 

e««(  el  +  eu  >/2  .0 

0024 

24  1 

cent inue 

0025 

24  i 

d3*en 

0024 

return 

0027 

end 
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0001 

0002 

0003 


UNIX 

0001 

0002 

0003 


UNIX 

0001 

0002 

0003 


UNIX 

0001 

0002 

0003 


UNIX 

0001 

0002 

9003 


real  function  d$< xd/n< ref / tref > 

46-0  .0 

end 
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real  function  d7< xd#n* ref / tref > 

d7«0.0  • 

end 
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real  function  d8( xd*n< ref • tref ) 

48-0.0 

end 
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real  function  d9( xd«n« ref < tref > 

49-0  .0 
end 
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real  function  dlO<xd« n »ref # tref > 

d 10-0 . 0 

end 
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DPE8N  .  FOR 

ft  library  contains  signal  waveforms  and  noise 
generator  functions  required  by  the  program 
DPEP 

for  all  signal  uave  form  functions  the  parameter 
js  represents  the  sample  order. 

For  all  noise  functions  the  parameter  an  (or  a) 
represents  the  scale  parameter  from  the  pregaran 
DPPRE . 
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0001 

root  function  s2ia; 

0002 

xa*l 00 . 0*a 

0003 

s2*xa 

0004 

roturn 

0005 

ond 

UNIX 

fortran  iv  v01-ll  source  listing 

pago 

00C1 

roal  function  s3<js) 

0002 

s3*0 .0 

0003 

and 
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0001  roal  function  s4(a> 

0002  s4*0.0 

0003  ond 
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0001 

0002 


roal  function  s5<a> 
s3«0 .0 
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0001  real  function  znl(an) 

c 

c  Nornal  Distribution 


c 


0002 

Integer  flag 

0003 

. 

if  <  f l ag  .eq.  12343)  go  to  420 

OOOS 

410 

xzn*2 . 0'  *an< 0/ 0>-  1.0 

0000 

yzn«2.0eran(0i0>-1.0 

0007 

zzn«xzn*xzneyzn*yzn 

0008 

if  <zzn  .ge.  1.0)  go  to  410 

0010 

zzn*sqrt<~2.0*alog<zzn)/zzn) 

0011 

znl8xzn*zzn*an 

0012 

yzn*yzn*zzn*an 

0013 

fUg-12343 

0014 

return 

0013 

420 

znl=yzn 

0018 

f lag*0 

0017 

return 

0018 

end 
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0001 


0002 

0003 

0004 


real  function  zn2(an) 
c 

c  Double  Exponential 
c 

zn2asign<l.*ran<0.0>-0.3>ean*-Q.707I0478aalog<ran<0 

return 

end 


0>  > 
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0001  real  function  zn3<an> 

c 
c 

c  A  Conponent  f ron  4  Lognornat  Distribution 
c 

0002  zn3*exp<znl< an  >>*cos< 6.263 1853* ran(0«  0> ) 

0003  return 

0004  end 
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0001 


0002 

0003 

0003 

0000 

0007 

0008 

0009 


real  function  zn4(an> 
c 

c  Nor na 1  Con tan mated  by  Nornat. 
c 

zz«ran<  0«  0  > 

if  <zz  .*#.  0 . 1 >  go  to  441 

zn4«znt<en> 

return 

441  zn4*znl<l  .0  > 
return 
end 
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0001 


reel  function  znS(en) 


c  Nor  no  l  Conton mated  by  Double  Exponential 
c 


0002 

n«rtn(  0/0) 

0003 

if  In  .go. 

OOOS 

zn3*zn2(en> 

0000 

return 

0007  431 

zn3*znl(l .0) 

0000 

return 

0009 

end 
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0001  real  function  zn6(a) 

0002  znt-C.O 

0003  tnd 
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0001  roll  function  zn?<a> 

0002  zn7«0.0 

0003  ond 
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0001  rtal  function  zn8(a> 

0002  zn8»0.0 

0003  ond 
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•  0001  roal  function  zn9(a> 

0002  zn9«0 . 0 

0003  ond 
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0001  roal  function  znlO<a> 

0002  znl0«0.0 

0003  ond 
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